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ABSTRACT 

Context. We study the jet kinematics of the blazar S5 0716+714 during its active state in 2003-2004 with multi-epoch VLBI obser- 
vations. 

Aims. We present a kinematic analysis of the large-scale (0-12 mas) jet of 0716+714, based on the results of six epochs of VLBA 
monitoring at 5 GHz. Additionally, we compare kinematic results obtained with two imaging methods based on different deconvolu- 
tion algorithms. 

Methods. The blazar 0716+714 has a diffuse large-scale jet, which is very faint compared with the bright compact core. Experiments 
with simulated data showed that the conventional data reduction procedure based on the CLEAN deconvolution algorithm does not 
perform well in restoring this type of structure. This might be the reason why previous kinematic studies of this source yielded am- 
biguous results. In order to obtain accurate kinematics of this source, we independently applied two imaging techniques to the raw 
data: the conventional method, based on difference mapping, which uses CLEAN deconvolution, and the generalized maximum en- 
tropy method (GMEM) realized in the VLB Imager package developed at the Pulkovo Observatory in Russia. 

Results. The results of both methods give us a consistent kinematic scenario: the large-scale jet of 0716+714 is diffuse and stationary. 
Differences between the inner (0-1 mas) and outer (1-12 mas) regions of the jet in brightness and velocity of the components could 
be explained by the bending of the jet, which causes the angle between the jet direction and the line of sight to change from ~ 5° to 

~ ir. 

Conclusions. We tested the performance of the two imaging methods on real data and found that they yield similar kinematic results, 
but determination of the jet component positions by the conventional method was less precise. The method based on the GMEM 
algorithm is suitable for kinematic studies. It is especially effective for dim diffuse sources with the average brightness of several mJy 
with bright point-like features. For the source 0716+714, both methods worked at the limit of their capability. 

Key words, galaxies: active - BE Lacertae objects: general - galaxies: BE Lacertae objects: individual: S5 0716+714 



1. Introduction 



The blazar 85 0716+714 at redshift z=0.31 (iNilsson et al.ll2008h 
is one of the most active sources of its class. It is highly 
variable on timescales from hours to months across all ob- 
served wavelengths. Intr a-Day Variability ( IDV) has been de- 
tected in this blazar (I Wagner et al.l 119961: IStalin et al.l 120061: 
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iMontagni et al.ll2006l: [Gupta et al.ll2009l) . T he IDV in optical an d 
radio bands are reported to be correlated (Wagner et al.l ll996h , 
which makes it the only candidate for a source with an intrin- 
sic origin of fast variation, which implies the presence of very 
small emitting regions within the source. At larger time scales 
this source has frequent flares in both radio and optical bands. 
Amp litudes of the flares i ncreased ab r uptly a fter the yea r 2003 
(e.g. iRaiteri et al.1 120031: iNesci et al.1 120051: iBach et al.1 12QQ7I: 
iFuhrmann et al.ll2Q08l and references therein). In 2007 and 2008 
it was also detected at very high-energy gamma rays with the 
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MAGIC telescope (iKranich & the MAGIC collaboratioiill2QQ9l : 
lAnderhub et aDl2QQ9^ 

At the milliarcsecond scale (VLBI observations) 0716+714 
has a core-dominated structure with a bright point-like core and 
a faint diffuse jet pointing at the position angle of ~ 20°. At 
the scale of arcseconds (VLA observations), the source struc- 
ture is also core-dominated, its jet is double-sided, located along 
the position angle of ~ 100° (misaligned with the r nilliarc second 
scale jet by ~ 80°), and has extended diffuse lobes (IWagner et alJ 
11996) . 

The long-term kinematics of this source have been exten- 
sively studied, and two mutually exclusive models have been 
suggested: one implies relatively fast outward motion of the su- 
perluminal components in the jet with speeds ranging from 5c 
to 16c ("fast scenario". iBach et aP 12005 ). and the other claims 
the jet to be stationary, with significant non-radial ("pr ecession") 
motion of the jet as a whole ("stationary scenario" Br itzen et"aD 
2009). Both these scenarios are based on the reanalysis of his- 
torical VLBI observations that span the years 1992-2001 ("fast") 
and 1992-2006 ("stationary") at frequencies from 5 to 43 GHz. 

Our analysis of the inner jet kinematics of 0716+714 at fre- 
quencies of 22, 43, and 86 GHz also revealed a fast-moving jet 
with appar ent speeds from 8 c to 19c on the spatial scale from 
to 1.5 mas (iRastorgueva et al.|[2009l) . Here we present results of 
our short-term kinematic study at 5 GHz. The data set covers six 
months, and the source was observed nearly every month. With 
this dense data sampling we are able to study the fast structural 
changes in the jet of 0716+714. 

The standard method of VLBI data reduction, which is used 
by the vast majority of researchers in the field, includes the 
Difmap package ( Shepherd 1997) as an imaging facility, which 
became conventional because of its simplicity and accessibil- 
ity. However, Difmap has its downsides, and the main disad- 
vantage of this package, w hic h uses the standar d CLEAN al- 
gorithm by Hogbom (^ 19741) (iShepherd et al.lll994[). is a pr oblem 
with the reconstruction of extended emission JComwe U|[l983). 
Tests performed on the simulated data have shown that the stan- 
dard CLEAN, and Difmap in general, are especially inaccu- 
rate in the reconstruction of the sources whose structure could 
be described as bright point-like knots embedded in a region 
with faint extended emission (e.g. iBaikoval [20071: iLawson et al.l 
|2004|) . The morphology of 0716+714 at lower frequencies (5- 
1.6 GHz) is consistent with this description. Therefore, a con- 
ventional approach to this source may have yielded erroneous 
results. In order to make sure that the structure of 0716+714 is 
correctly reconstructed, we independently applied both the gen- 
eralized maximum entropy method (GMEM) and the conven- 
tional CLEAN to the same data set. 

This paper is arranged as follows: section |2] describes the 
observations, section [3] contains a description of the a-priori 
data calibration and an introduction to the imaging methods. In 
section |4] we present and discuss the kinematic results. The fi- 
nal kinematic scenario, a comparison with previous studies, and 
possible explanations of the apparent speed differences are pre- 
sented in section O 

2. Observations 

We observed 0716+714 with the Very Long Baseline Array 
(VLBA) five times during the year 2004. Each observing ses- 
sion lasted nine hours and covered five frequencies (86, 43, 22, 
5 and 1.6 GHz) in dual polarization mode. Here we describe the 
results of the 5 GHz observations. The five epochs in our exper- 
iment were named in alphabetical order: A, Feb 10; B, May 03; 



C, Jun 18; D, Jul 29 and E, Aug 29, 2004. Additionally, we used 
one epoch from the unpublished polarization observation project 
of 37 northern BL Lac objects from the Kuehr & Schmidt sam- 
ple ( Kuehr& Schmidt) 1 1990) by D. Gabuzda and M. Mahmud, 
in our notation epoch M, Mar 22, 2004. Epoch M fills the gap 
between A and B of our experiment, making the time sampling 
almost uniform with an interval of approximately one month be- 
tween observations. 

3. Data reduction 

3.1. Calibration 

The data were correlated at the Socorro VLBA correlator. A- 
priori amplitude and phase calibration were made at Tuorla 
Observatory with NRAO's Astronomical Image Processing 
System (AIPS). A/P5"s standard VLBAPROC procedures were 
used for the amplitude and phase calibrations, scans on 3C279 
were used to calculate instrumental delay and phase residuals, 
and OJ287 was used as an instrumental polarization (D-term) 
calibrator. The subsequent procedure for our two methods is de- 
scribed below. 

3.2. Image reconstruction: the conventional method 

The first method we used was the conventiona l method, which 
is a combination of CLEAN deconvolution (iHogbomI 1 19741) 
with self-calibration and uv-plane Gaussian model fitting. This 
method is used by most of the astronomical community in stud- 
ies of outflow structure and kinematics. 

We used the Caltech Difmap package (IShepherdlll997l) for 
imaging, employing a CLEAN algorithm in a self-calibration 
loop: a model of the source, consisting of point-like clean- 
components is created with CLEAN, then the antenna gains are 
self-calibrated with the model, and the procedure is repeated. 
Self-calibrations were performed in the following way: several 
phase self-calibration loops were followed by an amplitude one, 
with the solution interval of amplitude self-calibration gradu- 
ally decreasing from the time of the whole observation to the 
data averaging time (0.5 min). Finally, we performed amplitude 
self-calibration with the solution interval equal to the averaging 
time. Then the final clean-component model was discarded, and 
a new source model was created. We used the Difmap proce- 
dure, Modelfit, to create a model consisting of circular compo- 
nents with a Gaussian distribution of flux, and fit it to the visibil- 
ity data. Except for one case, only circular Gaussian components 
were used to parameterize the data in order to reduce the number 
of free parameters. For all models, the brightest optically thick 
component at the southernmost end of the jet was considered 
to be the core, as in the previous research. Positions of other 
components were measured with regard to this. The resulting 
maps with the superimposed model components are presented in 
Fig. [T] and map parameters for all epochs are listed in Table [T] 
The errors of the position and the flux estimates of the compo - 
nents were calculated with the Difwrap package (lLovellll200Ql) . 
where the parameter to be analyzed (position, flux or FWHM) 
is varied around the best-fit value in a grid of a given step size 
and textitModelfit is run for each step. The result of this process 
is plotted as a distribution, and the user can visually inspect 
the model vs a uv-data plot or a residual map for each point of 
the grid. Visual inspection is based on the relative changes. 
Because the number of degrees of freedom of the model are 
not exactly known, this method of error estimate based on vi- 
sual inspection yields more reliable results than the standard 
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Fig. 1. Images of 0716+714 obtained using the conventional method: CLEAN maps, convolved with the averaged beam of 1.5x2.5 
mas, PA = -4.5°. Circles with crosses are the Modelfit Gaussian components superimposed on the image. The Gaussian model 
fitting was performed in the uv-plane. Axes are right ascension and declination in milliarcseconds, measured with regard to the 
core. 



Table 1. Map parameters for the conventional method. 



Epoch 



^ tot 

[Jy] 



Beam 
masx mas,[^] 



^ peak 

[Jy/beam] 



3o- 

[mJy/beam] 



2004.11 
2004.22 
2004.34 
2004.46 
2004.58 
2004.66 



1.46 
1.13 
1.06 
1.16 
0.90 
0.87 



2.5x1.5,-4.5 
2.5x1.5,-4.5 
2.5x1.5,-4.5 
2.5x1.5,-4.5 
2.5x1.5,-4.5 
2.5x1.5,-4.5 



1.41 
1.08 
1.01 
1.10 
0.84 
0.81 



0.33 
0.89 
0.36 
0.54 
0.72 
0.63 



analysis. In our case it was not possible to use the model vs a 
visibility plot, because we dealt with very faint components of 
the jet: changes in the model visibilities between diff'erent grid 
points were too small to be detected with the naked eye. Instead, 
we estimated the goodness of fit based on the appearance of the 



residual map and its changes from one gri d point to another . 
We applied an intuitive criterion proposed by ISavolainenI (l2006l) 
to determine the maximal acceptable deviation from the best- 
fit model: the diff'erence between a model and observed data is 
considered significant if the structure appearing in the residual 
map would look strong enough that it should be cleaned, if it 
had appeared during the imaging process. Because the flux of 
such a feature changes gradually from one grid point to another, 
it may be difficult to determine at which step the discrepancy be- 
tween the observed data and the model becomes large enough. 
In these cases we used an empirical '2(T-rule,' which we worked 
out by analysis of both the imaging process and the error esti- 
mate procedure: the flux is considered 'excessive' and should be 
'cleaned out' by placing a model component in its location if 
the rms noise in the area under scrutiny is twice as high as the 
average rms noise at the outskirts of the map. At the grid point 
where the noise level exceeds twice the rms, we assumed that the 
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parameter value had reached the limit. This deviation from the 
best-fit value was assumed to be the limiting error, A. A more de- 
tailed discussion of VLBI data error estimate with Difwrap and 
a compariso n with the method of variance about the best fit can 
be found in SavolainenI (|2QQ6|) . 

In most cases the Difwrap package yields asymmetric uncer- 
tainties for the model parameters, which leads to a bias in the 
resulting mo dels. This bias was corrected for with the method 
suggested by lD'Agostinil (l2QQ3l) . based on the probabilistic ap- 
proach to the bias correction: let us assume that the negative lim- 
iting error of the component's separation from the core r is A^^^r, 
and the positive limiting error is A^^^, and /^neg ^ ^pos- We cal- 
culated the unbiased limiting error as an arithmetic average of 
the two limiting errors: Aun=(^neg-^^pos)/'^- The standard devi- 
ation of r was calculated as crr=AunP- The unbiased value of 
the model parameter (in this case, separation from the core) is 
not equal to the result of the model fitting and is r^„=r-h(Apo^- 

^neg)/'^- 



3.4. Image reconstruction: the generalized maximum entropy 
method (GMEM) 

The eneralized maximum entropy method has been developed 
and implemented in the Pulkovo VLBI da ta reduction soft- 
ware package VLBImager by A. Baikova (lBaikovall2007h . The 
mapping technique realized in this package utilizes a self- 
calibration algorithm (Cornwell & Eomalont 19 9% combine d 
with the GMEM OFrieden & Bajkova 1994 iBaikoval l2007h . 
which is used as a deconvolution procedure. 

The GMEM is designed for the reconstruction of sign- 
variable functions, therefore it allows one to obtain unbiased so- 
lutions. The bias of the solution is one of the problems of the 
conventional MEM (Corn well et al.r i999), which may lead to 
a substantial n onlinear distort ion of the final image if the data 
contain errors (lBaykovalll995h . 

For the GMEM, the Shannon-entropy functional has the 
form 



3.3. Image reconstruction: the maximum entropy method 
(MEM), an attempt 

Another popular deconvolution algorithm in radio interf erome- 
try is the maximum entropy method invented by Jaynes (iJaynesI 
ll957), whi ch was firs t applied to image reconstruction by 
Frieden (Fr iedenlll972l) . The maximum entropy method is an 
extremely nonlinear method that leads to maximally smoothed 
images subject to constraints. The positivity of the sought-for 
solution forces super-resolution on bright isolated o bjects. 

There are many forms of entropy (Naraya n & Nityanandal 

Il986h . We used the A IPS task VTESS which uses the entropy 
functional in the form (IComwe lillT995h 



E(I) = - J I(x, y) ln[/(x, y)/M(x, y)]dxdy , 
I(x,y) > 0, 



(1) 



where M(x, y) is a default image that represents the optimum so- 
lution in the absence of data. Usually, the default image is taken 
to be flat or a lower resolution image. 

Given a form of entropy, the MEM defines the best image 
that maximizes E(I) subject to the data constraints: 1) fit the 
observed visibilities Vk to within the error term 6k'. F(I) = V)^+6y^; 
2) fit the total flux: J I(x, y)dxdy = Fq; 3) the rms error must be 

plausible: Z ^ = 

In AIPS, the Cornwell-Evans algorithm (IComwell & EvansI 

is used for the realization of the MEM. This algorithm 
uses a simple Newton-Raphson approach to optimize the func- 
tional ([B subject to the constraints upon the rms and total power 
enforced by the Lagrange multipliers. 

We made an attempt to use the MEM method to recover the 
faint emission of the 0716-h714 jet. We imported visibility data, 
self-calibrated in Difmap (see Section [T2b , to AIPS and used the 
task VTESS to create the final image. 

We used the task IMAGR to clean the central part of the im- 
age containing the bright core in order to avoid circular- shaped 
artifacts around the point-like source. Nevertheless, the image 
quality of this method was poor and did not allow us to correctly 
perform the kinematic analysis. Therefore, we conclude that the 
MEM method is not suitable for recovering the emission of the 
sources with bright point-like features and faint underlying dif- 
fuse emission. 



E{a) = - J (xP(t) ln(axP(t)) -h jc"(0 ln(ax"(0)) dt 
xP(t)>0, x"(0>0. 



(2) 



where xP(t) and x"^(t) are the positive and negative components 
of the sought-for image x(t), i.e. the equation x(t) = xP(t) - x"(t) 
holds, a > is a parameter responsible for the accuracy of the 
separation of the negative and positive components of the solu- 
tion x(t), and therefore critical for the resulting image fidelity. It 
is easy to show (lBaikovall2QQ7l) that solutions for xP(t) and x"(0 
obtained with the Lagrange optimization method are connected 
by the expression xP(t) • x^(t) = exp(-2 - 2 Ina^) = K(a), which 
depends only on the parameter a. This parameter is responsible 
for dividing the positive and negative parts of the solution: the 
larger a is, the more accurate is the discrimination. On the other 
hand, the value of a is constrained by computational limitations. 
The main constraint comes from the term in the optimized 
functional, which depends on the data errors. The larger a stan- 
dard deviation is, the higher value of a could be set. If data are 
very accurate, a lower value of a is needed. In practice, a is cho- 
sen empirically. In our case we had to compromise between data 
errors, which determine resolution of the resulting MEM solu- 
tion, and a need to divide the positive and negative parts of the 
solution as accurately as possible. It is fair to say that given fixed 
errors in the data, a maximum possible chosen value of a pro- 
vides us with the best possible resolution of the MEM-solution. 
In this work we used a = 100000. The algorithm takes into ac- 
count errors of the measurements. It implies a solution of the 
following two-dimensional problem of conditional optimization, 
given in the discrete form 



(3) 



m=l m 



N N 



k^l 1=1 



k=i 1=1 
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where k and / are the pixel numbers in a source map of the size 
N X N, m = l,...,Misa number of the visibility function mea- 
surement, a^^ and b^^ are the constant coefficients corresponding 
to the Fourier-transform, and are the real and imaginary 
parts of the complex visibility function (the Fourier spectrum of 
a source) and t/J^ and t/J^ are the real and imaginary parts of the 
additive visibility noise with normal distribution over zero value 
with known dispersion cr^n. 

As we can see from equation ©, the optimized functional 
has two parts: a Shannon-entropy functional and a functional 
that is an estimate of the difference between the reconstructed 
visibility spectrum and the measured data according to a cri- 
terion. This latter functional can be considered to be an addi- 
tional stabilizing term acting to provide a further regularization 
of the solution above what is possible with the entropy func- 
tional alone. Solution normalization, necessary for the MEM is 
obtained by limiting the total flux of the source: 



Numerical experiments show that the GMEM combined with 
the difl'erence mapping technique allows the restoration of the 
source structures consisting of bright compact features embed- 
ded in a weaker extended background with maximal accuracy 
(iBaikova 2007). 

After the a priori calibration, the raw visibility data were 
averaged in frequency (averaged IFs) in AIPS, then exported to 
Difmap, where they were rebinned, phase self-calibrated with 
a point source model with a 0.5 minute solution interval, and 
finally averaged in time (over a period of four minutes). After 
that the data was exported to VLBImager, where the automatic 
GMEM deconvolution - self-calibration loop iterations were 
performed until convergence. The size for the maps is chosen 
256 X 256 pixels with a pixel size of 0. 15 mas. The GMEM solu- 
tion convolved with a clean beam of ~ 2.4x 1 .5 mas, PA = -4.5° 
is presented in Fig. |2l and parameters for those maps are listed 
in Table O 



Table 2. Map parameters for the GMEM. The two tables corre- 
spond to two diff'erent sizes of the beam, with which the original 
GMEM solutions were convolved. 



Epoch 


S tot 


- Beam 


^ peak 


3cr 




[Jy] 


masx mas,['^] 


[Jy/beam] 


[mJy/beam] 


2004.11 


1.48 


2.4x1.5,-4.5 


1.42 


0.97 


2004.22 


1.08 


2.4x1.5,-4.5 


1.03 


0.96 


2004.34 


1.08 


2.4x1.5,-4.5 


1.00 


0.74 


2004.46 


1.15 


2.4x1.5,-4.5 


1.09 


0.99 


2004.58 


0.90 


2.4x1.5,-4.5 


0.82 


0.59 


2004.66 


0.86 


2.4x1.5,-4.5 


0.78 


0.67 


2004.11 


1.48 


0.5x0.5,0 


1.08 


0.33 


2004.22 


1.08 


0.5x0.5,0 


0.73 


0.29 


2004.34 


1.08 


0.5x0.5,0 


0.70 


0.24 


2004.46 


1.15 


0.5x0.5,0 


0.74 


0.30 


2004.58 


0.90 


0.5x0.5,0 


0.53 


0.19 


2004.66 


0.86 


0.5x0.5,0 


0.48 


0.19 



In the appendix we present two additional sets of maps to 
demonstrate the capabilities of this deconvolution method: the 
original GMEM solutions (Fig.O and the GMEM solutions con- 
volved with a circular clean beam of the size 0.5x0.5 mas, which 



were used for the modeling of the source structure (Fig. |7). The 
model-fitting of the maps was made in the image plane. The 
source structure was modeled using Gaussian elliptical compo- 
nents. Each of the components was described by the following 
six parameters: r, distance from the core (the core was deter- 
mined like in the conventional method and was assumed to be 
stationary); 0, position angle measured counter-clockwise be- 
tween the ordinate axis and the radius-vector: and 

^ mill 9 Sizes 

of the major and minor axes of a Gaussian ellipsoid at half maxi- 
mum amplitude and 6, the angle describing the orientation of the 
component ellipse, an angle between the ordinate axis and the 
major axis of the ellipse, measured counter-clockwise. Cross- 
identification of the components was based on the assumption 
that the core separation, the position angle and the flux of the 
component change smoothly from epoch to epoch. In order to 
obtain a model of the inner jet (r<1.5 mas) with higher accu- 
racy, we interpolated the brightness distribution to a finer grid: 
the pixel size in thit case was 0.0375 mas, which was four times 
smaller than that for the large-scale jet. The errors of the com- 
ponent core separation determination were estimated by the for- 
mula (.Fomalontl999.) 



Ar : 



0-0 



2^ 



peak 



where cr is the rms after the subtraction of the model from the 
source map, 5'peak is the peak flux of the component, is the size 
of the component at the amplitude half maximum (in our case we 
assume = amaj). We point out that that this formula likely un- 
derestimates position errors, which is clear from the comparison 
with the Difwrap errors (See Fig. |3). 



4. Results 

Our aim was to study the structure and kinematics of the large- 
scale (from 1 to 12 mas) jet, and compare it with the re- 
sults for the inner (< 1 mas) jet, which were published in 
iRasto rgueva et al.' ('2009'). 

We applied a conventional method and a GMEM to the same 
data set to compare their kinematic results. For each of them, 
we calculated the core separation of the components as a func- 
tion of time, the proper motion yu, and the apparent speeds jSapp 
according to the formula 



/^app — 



c{l+z) ' 



where the luminosity distance Jl is calculated with the analytical 
formula derived by Pen (1999^) for a flat cosni o logy. The same 
analytical expression was used by iBach et aP (l2005h for their 
kinematic analysis of this source, therefore our results could 
be compared directly without conversion. For this cosmological 
model, a proper motion of 1 mas/yr is equivalent to the appar- 
ent speed of 19.26c. The flat Universe model with the Hubble 
constant of //q = 71km/s/M pc and with = 0.3 and Qa - Q-7 
was used. We are aware that Wickramasin she & Ukwattal (l2010l) 
recently obtained a new and more precise analytical expres sion 
for the luminosity distance, but for z=0.3 the difference with lPenI 
(1999) is negligible (1 mas/yr corresponds to the apparent speed 
of 19.22c). 

The results of the Gaussian model-fitting and component 
identification are presented in Table |5] (conventional method) 
and in the Table [6] (GMEM). The original GMEM- solutions are 
shown in Fig.[6l 
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Fig. 2. GMEM solutions for each epoch, convolved with the averaged CLEAN beam: 1.5x2.5 mas, PA = -4.5°. Gray ellipses 
are the Gaussian components superimposed on the image. The Gaussian model-fitting was performed in the image plane for the 
original GMEM solutions (Fig.O convolved with the beam of 0.5 x 0.5 mas (Fig. [7]). Axes are right ascension vs declination in 
milliarcseconds, measured with regard to the core. 



4.1. Component cross-identification between different 
methods. 

4.1 .1 . Positions and apparent speeds of components 

Based on core separation and fluxes, we identified four nearly 
stationary components in the large-scale jet of 0716-h714: one 
at a separation of about -1.5 mas (C3 - A3), one at -3.0 mas 
(C4 - A4), one at -4.6 mas (C5 - A5) and one at -12 mas (C6 - 
A6). Their proper motions and apparent speeds can be found in 
Tables [3] and m A comparison of their positions, obtained with 
difl'erent methods, and linear fits are shown in Fig. [3] The com- 
ponents C1-C6, found by the conventional method, have a larger 



scatter around the mean position and larger positional errors 
than those of the GMEM components (compare Tables [3] and|4l), 
which complicated the component cross-identification and made 
a choice of the kinematic scenario ambiguous. However, the 
GMEM method confirmed the stationary scenario. 

However, the component CI -2, determined by the conven- 
tional method, and the component Al found by the GMEM in 
the inner jet have a clear outward motion: - 10c and - 2c, re- 
spectively (Tables [3] and |4]). The diff'erence in apparent speed 
could be explained by the fact that the GMEM determined an ex- 
tra stationary component at the separation of about 1 mas from 
the core (A2), while the conventional method did not. On the 
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Table 3. Kinematic scenario obtained with the conventional 
method. The table contains the component name (enumeration 
starts from the vicinity of the core), the average distance from 
the core in mas, the component proper motion in mas/yr and ap- 
parent speed in the units of speed of light. 



Comp. 


Core dist. 




Papp 


name 


[mas] 


[mas/yr] 


[c] 


Cl-2 


0.64+0.14 


0.55+0.11 


10.6+2.2 


C3 


1.9+0.4 


0.62+0.41 


11.8+8.0 


C4 


3.2+0.4 


0.44+1.4 


8.6+26.5 


C5 


5.7+0.9 


-1.37+2.7 


-26.4+52.2 


C6 


11.6+0.6 


-0.27+1.2 


5.2+23.6 



Table 4. Kinematic scenario obtained with the GIMEIM. The ta- 
ble contains the component name (enumeration starts from the 
vicinity of the core), the average distance from the core in mas, 
the component proper motion in mas/yr and apparent speed in 
the units of speed of light. 



Comp. 


Core dist. 




Papp 


name 


[mas] 


[mas/yr] 


[c] 


Al 


0.49+0.02 


0.12+0.05 


2.4 +0.9 


A2 


0.93+0.06 


-0.24+0.40 


-4.6+7.1 


A3 


1.48+0.13 


-0.64+0.03 


-12.4+0.6 


A4 


3.3+0.2 


0.13+0.32 


2.4 +6.2 


A5 


4.7+0.4 


-1.37+0.37 


-26+7 


A6* 


11.9+0.8 


-0.45+1.13 


-9+20 



* fit performed without taking into account errors of the component 
position estimate 
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Fig. 3. All components obtained with the conventional and 
GIMEIM methods are shown at the same plot. 



other hand, A2 is only found at the last four epochs, and is sep- 
arated from the component Al by only -0.5 mas, and the com- 
ponent Cl-2, determined by the conventional method, is located 
between them (see Fig.[3l). Therefore, we believe that the compo- 
nent Cl-2 is resolved by GIMEIM in the four last epochs, making 
up two components Al and A2, and this effect causes the differ- 
ence in the apparent speed. 

Note that the kinematic s of the outer jet are consistent with 
the results of iBritzen et"aD (12009.) . 

4.1.2. Core flux 

In this section we compare the ability of the two methods to re- 
store the flux of the source. In the optically thin jet, the fluxes 
of the components found by the conventional method and by 
the GIMEIM agree within two sigma. At the same time, the per- 
formance of the two methods in the core region is significantly 
different: the conventional method places one very bright com- 
ponent Cl-2 (the flux density is about 0.1 Jy, so it is bright com- 
pared to other jet components) next to the core, while the GIMEIM 
method finds two components, Al and A2, with the flux density 
comparable to the other components in the jet (on the order of 
0.01 Jy). However, if one adds up the flux densities of the core 
and nearby components (CO+Cl-2 and A0+A1+A2), the results 
of the two methods would agree within two sigma for all six 
epochs. We attribute this diflTerence to the fact that the GIMEIM 
restores the bright point-like core accurately, whereas the con- 
ventional method, because of significantly lower resolution at 5 
GHz, is likely to mix up emission of the core and of the nearby 
components. 



5. Discussion 

5.1. Inner jet kinematics: comparison of the 22 and 5 GHz 
data 

The inner jet of 0716+714 is well resolved at higher (> 5 GHz) 
frequencies and reveals a clear outward motion, although the 
emission fades away at ab out 1 mas as the spectral peak moves 
to lower frequencies (see iRastorgueva et al T l2009h . At 5 GHz, 
only one or two components are resolved in the inner jet region, 
depending on the method, and they also have a clear outward 
motion. Because the data at higher frequencies were reduced by 
the conventional method, we discussin in this section only the 
component Cl-2. Its apparent speed can be found in Table [3] 
In Fig.|4]we show the motion of the component Cl-2 at 5 GHz 
together with components C5 and C6 at 22 GHz. They are obvi- 
ously located in the same area in the inner jet and have speeds of 
a similar values: the apparent speeds of the components C5(22 
GHz) and C6(22 GHz) are 21.3+2.5 and 18.9+0.8 respectively, 
which agrees with the speed of CI -2(5 GHz) within 3cr. The 
diflTerence in the apparent speeds of the inner jet components, 
detected at diflTerent frequencies, is likely caused by a resolution 
eflTect. 

5.2. Jet bending and beaming 

The jet components that were found beyond 1 mas from the core 
at both 5 and 22 GHz have a very high dispersion in the appar- 
ent speeds, and in most cases move slower than the inner ones 
(see Table O. Therefore, we observe two kinematically diflTer- 
ent regions in the jet of 0716+714: a fast moving inner jet (< 1 
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Fig. 4. 22 GHz inner jet components C5 and C6 from 
iRas torguev a et alJ (l2009l) and the component CI -2 at 5 GHz co- 
incide spatially and have similar motion. 



mas) and a slowly moving outer jet with nearly stationary com- 
ponents (> 1 mas). The orientation of the inner and outer parts 
of the jet on the sky plane are also different: the jet changes its 
position angle from ~ 25° (inner region) to ~ 15° (outer region). 
The bend occurs at the distance of about 1 mas from the core, 
see Fig. [51 In addition to that, we found that the flux of the jet 
components in the inner jet is systematically higher than the flux 
of the outer jet components. Based on this we propose that those 
efl'ects could be caused by the jet bending away from the line of 
sight at the apparent distance of 1 mas from the core. We tried 
to find out whether it is possible to explain the observed eff'ects 
by the jet bending, without additional assumptions of jet param- 
eters. If the jet bends away from the line of sight, the beaming 
factor 6 decreases, which causes a decrease in the apparent lumi- 
nos ity of the jet, beca use it is proportional to 6^^^ for a smooth 
jet (ICohen et al.ll2QQ7l) . where a is a spectral index {S ~ v^). The 
apparent speed of the jet also changes with the change of the an- 
gle to the line of sight, for 0716-h714 it decreases. We estimated 
how much the flux of the jet components change from the inner 
jet to the outer, and compared these results with the theoretical 
predictions for jet bending. We assumed the jet to be a smooth 
collimated flow with a constant plasma speed (and therefore con- 
stant Lorentz factor F). This assumption is indirectly supported 
by the fact that the jet components are well aligned (see Fig.O, 
and the structure of the jet on VLBI maps is smooth and con- 
tinuous. We used the results of the conventional method for this 
analysis. 

For each epoch we calculated the ratio of the average com- 
ponent flux at 5 GHz in the outer part of the jet to the average 
flux at 22 GHz of the inner ones: rj =< lo > / < It >. Before cal- 
culating the flux ratio r/ for each epoch, we corrected the spec- 
tral slope of the 5 GHz fluxes in order to compare them directly 
with the 22 GHz data. The 5-22 GHz spectral index a was es- 
timated for the inner jet components CI -2 (5 GHz) and C5 and 
C6 (22 GHz) for each epoch separately. Individual values of the 
flux ratios were averaged, which gave < rj >= 0.11 ± 0.08. 
Assuming that the flux change is caused by jet bending, we 
can calculate the ratio of the beaming factors of the jet after 
and before the bent using the flux ratio we have just obtained: 
< rs >= = 0.24 ± 0.12. These values were also calcu- 
lated separately for each epoch, using the individual values of the 
spectral index, and then averaged. Errors of the averaged value 
consist of the standard deviation of averaging the ratios between 
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Fig. 5. Components of the inner and outer jet on the sky plane. 
The change of the orientation occurs at the distance of about 
1 mas from the core. The components of the jet are well aligned, 
which is an indication of the collimated jet. 

epochs and the standard deviation of averaging the component 
fluxes for each epoch, added in quadrature. 

The apparent kinematics of the relativistic blazar jet are 
determined by four parameters: the apparent speed /3app, the 
Lorentz factor F, the beaming factor S and the angle of the jet 
to the line of sight Q, any two of which could be used to find two 
others fe.g. lCohen et al.ll2007l) . In the inner and outer parts of the 
jet, the Lorentz factor F is the same, but the other three param- 
eters diff'er: ©2 > ©i, and therefore Papn 2 < ffappi and 62 < • 
From the 22 and 43 GHz data analysis (Rastorgu eva et al.|[2009h 
we estimated the kinematic parameters for the inner jet: the max- 
imum apparent speed is Pappy^ax ~ ^^-^ corresponding 
minimum Lorentz factor is Ymin = 19.5. The corresponding min- 
imum angle of the jet to the line of sight is &mins = 3.0°. In the 
outer jet, the maximum apparent speed was jSapp - 10. Taking 
into account the above mentioned observed quantities and as- 
suming that the jet has a constant speed, we can calculate its 
angle to the line of sight before (0i) and after (©2) the bending 
as well as other kinematic parameters of the jet. 

We can consider the beaming factor of the jet as a function 
of the initial (inner jet) angle to the line of sight ©i, which is 
larger than the &min = 3.0°, but smaller than the maximal value 
of &max = 3.rcsm(2/3app/(l + Plpp) = 5.9° (aU possible values lie 
in the lower part of a /3app - S diagram). The ratio of the beaming 
factors rs is a monotonic function of ©i and changes from 0. 15 to 
0.26, while ©i changes within its limits, and reaches the value of 
0.24 at ©1 = 5.2°. Therefore, the observed change of kinematics 
is possible if the Lorentz factor of the jet is F ^ 30.0, and at 
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the apparent distance of 1 mas from the core the jet changes its 
direction, with the angle to to the line of sight increasing from 
~ 5° to ~ 1 1°. At the same time the beaming factor is changing 
from ~ 7 to ~ 2, respectively. 

The errors of the observed flux factor ratio determination, 
which give us a beaming factor ratio, are fairly large, and precise 
calculations of the angle to the line of sight and beaming factor 
are not possible. The values we obtained are only estimates, and 
we can only say that the observed change of the jet apparent 
speed and brightness is possibly caused by jet bending. Also, the 
proposed beami ng factor for the inne r jet is consistent with the 
value found by iHovatta et al.l (l2009l) from the radio-frequency 
variability. 



6. Conclusions 

We studied the large-scale jet kinematics of the blazar S5 
0716+714 by analyzing six epochs of VLBA observations at 5 
GHz. Two imaging methods were applied to this data set: a con- 
ventional method (CLEAN self-calibration loop in Difmap com- 
bined with Modelfit) and GMEM (Pulkovo VLBImager package 
combined with the image plane model-fitting). 

6.1. Performance of the methods 

S5 0716+714 has a bright point-like core and a faint diff'use jet. 
The emission fades down very quickly in the vicinity of the core, 
and the CLEAN algorithm, producing the model that is com- 
posed of the discrete point sources, is not adequately restoring 
this emission region. In the Gaussian model, some fraction of 
the core emission was added to the inner jet component, and the 
kinematic results were also ambiguous because of the large un- 
certainties of the positions of jet components. 

The GMEM method delivered a smooth image of the source, 
which restored both the bright core and the diff'use jet. However, 
the jet of 0716+714 is very faint, therefore a convolution of the 
original solutions with the small (0.5x0.5 mas) beam was needed 
in order to obtain accurate Gaussian models of the jet structure. 
The original GMEM solutions convolved with the CLEAN beam 
in general resemble the CLEAN images (Fig.|2]and[T|). 

We recommend GMEM in combination with the diff'erence 
mapping technique for restoring the structure of compact AGN 
and other objects that have a diff'use structure with bright point- 
like features. In addition, the original GMEM solutions have a 
higher resolution than other methods, and in principle they could 
be used directly to derive kinematics of th e source, provided that 
the je t features are bright enough (e.g., iBajkova & PushkarevI 
I2008h . 

6.2. Structure and kinematics of the jet 

The large-scale jet of 0716+714 is faint, diff'use and smooth, 
without prominent brightness enhancements. 

Both the conventional method and GMEM deliver similar 
kinematic results: the large-scale jet (1-12 mas) is mostly sta- 
tionary, with a large scatter of the component positions. 

The position angle of the jet on the sky plane changes by 
~ 10° at about 1 mas from the core. 

Large-scale jet components are systematically slower and 
fainter than those of the inner (<1 mas) jet components. 

We attribute these last eff'ects to the bending of the jet at the 
apparent distance ~ 1 mas from the core from ~5°to~ll°. 
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Fig. 6. Original GMEM solutions. 
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Fig. 7. GMEM solutions, convolved with the beam of size 0.5x0.5 mas. 
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Table 5. Results from the CLEAN and Gaussian model fitting and component parameters. S peak is the peak flux density, r and PA 
are the distance and position angle measured from the core and agauss is the EWHM of the circular Gaussian component. ID is a 
component identification in one of the three possible scenarios, see Section|4l 
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Table 6. Results from the GMEM method imaging and model fitting with component parameters. S peak is the peak flux density, r 
and PA are the distance and position angle measured from the core, Umajor and aminor are the major and minor axes of the elliptical 
Gaussian components and ^ is a position angle of the ellipse. ID is a component identification in the stable jet scenario, see SectionlH 



Epoch 


Flux / 


T 


PA 


^ma jor 


^minor 


e 


ID 




rmJvl 


[mas] 


L J 


[mas] 


[mas] 


L J 




2004. 1 1 


1414.5 








0.61 


0.54 


8.99 


Core 




7.4 


486+0 006 


28.56 


0.61 


0.51 


138.14 


Al 




12.1 


1 664+0 007 


31.18 


0.75 


0.66 


-30.33 


A3 




3.0 


2 502+0 01 3 


5.24 


1.22 


0.76 


24.61 


X 




6.1 


3 470+0 043 


13.44 


1.20 


0.83 


117.53 


A4 




7.0 


5 3Q8+0 342 


12.72 


2.37 


1.94 


279.91 


A5 




1.1 


12 542 +0 679 


15.29 


5.57 


1.39 


53.18 


A6 


2004 22 


1011.8 








0.65 


0.55 


203.02 


Core 




19.9 


486+0 002 


22.08 


0.64 


0.57 


33.93 


Al 




14.9 


1 ^33+0 014 


25.75 


1.10 


0.89 


45.55 


A3 




8.9 


2 984+0 060 


14.01 


1.71 


0.92 


98.89 


A4 




3.8 


4 (S35+0 140 


15.09 


1.71 


0.90 


103.87 


A5 




3.5 


11 156 + 021 


8.65 


0.99 


0.62 


24.58 


A6 


2004.34 


945.6 








0.64 


0.55 


21.41 


Core 




28.7 


462+0 001 


26.63 


0.59 


0.52 


-60.54 


Al 




19.4 


986+0 004 


28.59 


0.87 


0.60 


34.05 


A2 




8.5 


1 483+0 016 


24.74 


1.15 


0.66 


145.82 


A3 




8.5 


3 315+0 022 


13.14 


1.14 


0.92 


-202.12 


A4 




3.2 


4 775+0 110 

1 • / / ^ _!_ Vy • 1 1 Vy 


19.86 


1.69 


0.77 


41.74 


A5 




3.2 


7 532+0 292 


19.04 


2.85 


0.71 


25.94 


X 




5.8 


13 10 + 666 

L ^) , 1 Vy _!_ vy.v/vyvy 


18.03 


4.99 


1.31 


38.75 


A6 


2004.46 


961.0 








0.65 


0.55 


21.80 


Core 




29.1 


489+0 001 

vy • 1 yJ y _ !_ Vy . Vy Vy 1 


26.83 


0.56 


0.53 


112.24 


Al 




32.5 


916+0 003 

Vy . V 1 Vy _!_ Vy . Vy Vy 


27.80 


0.96 


0.61 


-144.55 


A2 




6.6 


1 549+0 038 

\^ ,mJ^ y — !— Vy . Vy Vy 


16.70 


1.39 


0.69 


134.55 


A3 




6.1 


3 399+0 042 

^ m ^ y y — !— Vy . Vy i 


14.19 


1.30 


0.81 


-170.79 


A4 




3.5 


4 601 +0 240 

Vy Vy 1 _!_ Vy . .i— ^ i Vy 


20.80 


2.21 


0.91 


118.47 


A5 




11.8 


11.11 1.711 


16.25 


8.39 


2.36 


159.12 


A6 


2004 "SS 


711.7 








0.66 


0.56 


24.10 


Core 




36.0 


489+0 001 

Vy • 1 Vy y _ !_ Vy . Vy Vy 1 


25.70 


0.57 


0.56 


90.70 


Al 




38.5 


978+0 001 


27.75 


0.86 


0.60 


28.13 


A2 




6.9 


1 360+0 009 

1 . ^ Vy Vy _!_ Vy . Vy Vy V 


17.91 


0.94 


0.59 


6.29 


A3 




7.1 


3 367+0 025 


15.00 


1.29 


0.88 


-25.33 


A4 




4.5 


4.527+0.141 


22.39 


2.22 


1.09 


115.54 


A5 




13.3 


11.33 + 2.216 


18.12 


8.45 


3.61 


203.37 


A6 


2004.66 


712.2 








0.70 


0.56 


19.46 


Core 




27.4 


0.523+0.001 


21.35 


0.60 


0.55 


98.63 


Al 




7.9 


0.852+0.003 


21.66 


0.57 


0.53 


38.70 


A2 




25.8 


1.300+0.003 


25.19 


0.91 


0.66 


-150.48 


A3 




4.8 


3.558+0.108 


16.51 


2.07 


1.02 


-52.99 


A4 




4.8 


4.189+0.083 


15.81 


1.92 


0.89 


33.07 


A5 




1.4 


6.886+0.062 


19.18 


1.06 


0.62 


196.45 


X 




9.3 


12.40+0.794 


12.37 


6.83 


1.31 


59.65 


A6 



